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DM particles are the inert Higgs boson, % and (j). It is shown how the allowed parameter space, 
obtained previously in the absence of x an d (j), changes. The semi-annihilation process in this model 
produces monochromatic neutrinos. The observation rate of these monochromatic neutrinos from 
the Sun at IceCube is estimated. Observations of high energy monochromatic neutrinos from the 
Sun may indicate a multi-component DM system. 
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I. INTRODUCTION 



Recent astrophysical observations [lH3| have made it clear that most of the energy of the 
universe consists of dark energy and cold dark matter (DM), and their portions are very 
well fixed by these observations. While the origin of dark energy might be the cosmological 
constant of Einstein, the origin of cold DM cannot be found within the framework of the 
standard model (SM) of elementary particles. Moreover, we do not know very much about 
the detailed features of DM at present, even if the origin of DM should be elementary 
particles. Currently, many experiments are undertaken or planned, and it is widely believed 
;hat the existence of DM will be independently confirmed in the near future (see, for instance, 

1ft. 

A particle DM candidate can be made stable by an unbroken symmetry. The simplest 
possibility of such a symmetry is a parity, Z 2 . Whatever the origin of the Z 2 is, the lightest 
Z 2 odd particle can be a DM candidate if it is a neutral, weakly interacting and massive 
particle (WIMP) (see Q for a review). There is a variety of origins of the Z 2 . R parity 
in the minimal supersymmetric standard model (MSSM), which is introduced to forbid fast 
proton decay, is a well-known example (see 0] for a review). In this paper we consider a 
universe consisting of stable multi DM particles [zj]- 28]. A multi-component DM system can 
be realized if the DM stabilizing symmetry is larger than Z 2 : Zn (N > 4) or a product of 
two or more Z 2 s can yield a multi-component DM system. 1 In a supersymmetric extension 



of the radiative seesaw model of [3l| . for instance, a. Z 2 x ZL symmetry appears, providing 
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various concrete models of multi-component DM systems 

In a multi-component DM system, there can be various DM annihilation processes that 
are different from the standard DM annihilation process 32]- 37|, DM DM — > XX, where 
X is a generic SM particle in thermal equilibrium. Even in one-component DM systems, 
the non-standard annihilation process, the co-annihilation of DM and a nearly degenerate 



unstable particle 38], can play a crucial role in the MSSM 39]. The importance of non 



standard annihilation processes such as DM conversion 
DM 
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and semi-annihilation of 



22l | in two-component DM systems for the temperature evolution of the number 
density of DM has been recently reported. 

If {Z 2 ) e is unbroken, there can exist at least K = I stable DM particles. In a kinematically 
fortunate situation, 2 e — 1 stable DM particles can exist; for i = 2 there can be maximally 
K — 3 stable DM particles. Any one-component DM model can easily be extended to a 
multi-component DM system. The allowed parameter space of a one-component DM model 
can considerably change, as it has been recently found in Q (see also 0]), even using 
a crude approximation of a DM conversion process in a supersymmetric extension of the 
radiative seesaw model. 

In Sect. II, after we will have outlined a derivation of the coupled Boltzmann equations 



1 .Z3 allows only one-component DM systems. Refs. 



29 



discuss models with Z% 



2 



that are appropriate for our purpose, we will consider Active two- and three-component DM 
systems and analyze the effects of non-standard annihilation processes of DM. We will also 
address the question, within the framework of thermally produced DM (32! , whether one 
can obtain a large boost factor in a three-component DM system that is needed to explain 
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the positron excess observed by recent experiments 41J- 44J . We will find that in a 
fortunate situation a boost factor of O(10) can be obtained, which can be combined with 
the known mechanisms j^sj-js^], to obtain a desired enhancement effect. 
In Sect. Ill we will extend the radiative seesaw model of 



31 



by adding an extra Major ana 
fermion x an d an extra real scalar boson so as to obtain Z 2 DM stabilizing 

symmetry. Apart from the presence of 0, the Higgs sector is identical to that of 5lM 53 ] . 
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This model shall show how the allowed parameter space, which is obtained in 
the assumption that the lightest inert Higgs boson is DM, can change. Indirect detection of 
DM, in particular of neutrinos from the annihilation of the captured DM in the Sun 541]- 



63| are also discussed. We will solve the coupled evolution equations of the DM numbers 



in the Sun, which describe approaching equilibrium between the capture and annihilation 
(including conversion and semi-annihilation) rates of DM, and estimate the observation rates 
of neutrinos. Due to semi-annihilations of DM, monochromatic neutrinos are radiated from 
the Sun. Our conclusions are given in Sect. IV. 



II. BOLTZMANN EQUATION AND TWO- AND THREE-COMPONENT DM 
SYSTEMS 

A. Boltzmann equation 

Here we would like to outline a derivation of the Boltzmann equation that we are going to 
apply in the next section. We will do it for completeness, although the following discussion 



partially parallels that of |16|] (see also |2lJ). We start by assuming the existence of K stable 
DM particles x% with mass mj. None of the DM particles have the same quantum number 
with respect to the DM stabilizing symmetry. All the other particles are supposed to be in 
thermal equilibrium. Then we restrict ourselves to three types of processes which enter the 
Boltzmann equation: 

Xi X, Xi X[ , (1) 
Xi Xi Xj Xj (DM conversion) , (2) 
Xi Xj ^ Xk Xij k (DM semi-annihilation) , (3) 

where the extension to include coannihilations and annihilation processes like x% + Xj ^ 
Xk + Xi is straight forward. See Fig. [TJ for a depiction of DM conversion and DM semi- 
annihilation. We denote the phase space density of x% by fi(E{, t) and its number density by 
riiit) = ((yf/(27r) 3 ) / d 3 pifi(Ei,t), where g stands for the internal degrees of freedom. Then 
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FIG. 1: Dark matter conversion (left) and semi-annihilation (right). 



the density satisfies the Boltzmann equation (see, e.g., 35|), which we will not spell out 
here. Instead, we make the replacement 

t = 0.301g; 1/2 M PL T- 2 (4) 

during the radiation-dominated epoch, where t is the time of the comoving frame, g* is the 
total number of effective degrees of freedom, T and Mpl are the temperature and the Planck 
mass, respectively. Further, we use the approximation 

fi(Ej,t) Uiify 

fi(Ei,t) n t {t) ' 1 ] 

where fi(Ei,t) ~ exp(— Ei/T) and fii = (g/(2n) 3 ) J d 3 pifi(Ei,t) are the values in equilib- 
rium, and we ignore the chemical potential. Then the collision terms in the Boltzmann 
equation can be written as 

-(PSl)\M(ii;XiXi)\ 2 ^- (n ini -mi) 

riiUi 

-E(PSI) |M(iz;ii)| 2 # ( n,ii, - ^11,11) 



i>j 



1 2 fjfj ( - - 

ii ,ii : V n,n, • ' 



+ V(PSI) \M(jj:n)\--^-^ l rijTij n.,11 



-£(PSI) \M(ij;kX ijk )\ 2 ^- (n t n, - —n,n 3 

+ £(PSI) \M(jk;iX jkl )\ 2 ^ (n 3 n k - -n,n fc ) , (6) 
~y rijUk V rii ) 

where PSI stands for "phase space integral of (2tt) 4 <5 4 (momenta) x" , M is the matrix element 
of the corresponding process, and we have assumed that 

rrii > rrij for i > j and m^, mx'., m x l]k « mi for all i, j, k, I . (7) 
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Using the notion of the thermally-averaged cross section 

1 



< a(ii;XiXi)v > = —PSIlMfaXiX^fifi , (8) 



n,n. 



and the dimensionless inverse temperature x = fi/T, we obtain for the number per comoving 
volume Yi = m/s: 



( ' Y ' -0.264 g\l 2 [^p] { <a(U;X l Xl)v> (Y^-Y^ 



X' 

YY \ w , / yy _ _ 



dx 

■ X TO - £ttTO)- X < ( r(ii;n>>(y i y i - ^y^- 

+]T <t7(ij; fcx - fc ) v >(y i y j - ^yy,) -J2 ix ife >>(y i y fc - &y fc ) } , (9) 

j,fe ^ Yk ' j,k ^ Yi ' 

where = (X^ «\ rl ) is the reduced mass of the system. To arrive at Eq. we have used: 
s = (2it 2 /45)g*T s , H = 1.66 x g x J 2 T 2 / 'M PL , where s is the entropy density, and H is the 
Hubble constant. 

We will integrate this system of coupled non-linear differential equations numerically. 
Before we apply the Boltzmann equation Q to a concrete DM model, we discuss below 
the cases of K = 2 and 3 simply assuming fictitious values of the thermally-averaged cross 
sections and DM masses rrii. 



B. Two-component DM system 



Before we come to one of our main interests, a three-component DM system, we first 
consider the K = 2 case with a Z 2 x Z 2 symmetry. In this case, there are three different 
thermally-averaged cross sections. No semi- annihilation is allowed due to Z 2 x Z 2 . 2 We 
further assume that there are only s-wave contributions to < av > and X{ (i = 1, 2) are 
massless while mi > m 2 : 



Eq. 



<a(U;X 1 X[)v> = <t 0) i x 1(T 9 GeV~ 2 , < a(22; X 2 X' 2 )v >= a 0>2 x 10~ 9 GeV~ 
<a(ll;22)t»> = a 0>12 x KT 9 GeV~ 2 . 

© then becomes 



(10) 



dY 1 
dx 



-0.264 g x J 2 



fiMi 



PL 



Xr 



{<a(ii;x 1 x[)v> (yy - yy 



i <vr ( n:22)r>(yy-^y 1 y 1 

J- 2 J- 2 



(11) 



In Refs. , the Z4 case is discussed in detail. In this case there exist two independent DM particles, 

because due to CP invariance the anti-particle is not an independent degree of freedom in the Boltzmann 
equation. Semi-annihilation is allowed in this case. 
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FIG. 2: Left: The relic abundance £l Xl h 2 (x) (black) and il X2 h 2 (x) (blue) as a function of x = fi/T = 
[(m^ 1 +rrq 1 )T]- 1 with ct ,i = 0.1 , a , 2 = 6 , <r ,i2 = 4.4(0) (solid (dashed)) , mi = 200 GcV, m 2 = 160 
GeV and g* =90. Right: The total relic abundance flxh as a function of co,i2 which parametrizes the size 
of the conversion xiXi ~^ XiXi- 



dYo 



dx 



-0.264 gi 



1/2 



PL 



{<0(22;X 2 X 2 )*;> (r 2 F 2 - F 2 F 2 



<a(ll;22)v>[Y 1 Y 1 



Y 2 Y 2 



Y 1 Y 1 



(12) 



We consider the case in which the size of the DM conversion and the standard annihila- 



tion are of similar order (see also 2l|). In Fig. |2] (left) we show the evolution of the fraction 



tively, where we have used O ,i = 0.1 , cr ,2 = 6 . 
GeV, m 2 = 160 GeV and g* = 90, and x = \ijT 



of critical densities, Q xi h 2 (x) (black) and Q X2 h 2 (x) (blue), contributed by Xi an d X2, respec- 

a 012 = 4.4(0) (solid (dashed)) , m 1 = 200 
= [(^If 1 + m 2 )T\~ 1 . As we see from Fig. [2] 
(left), at o"o,i 2 = (i.e., no DM conversion (j2J)), the density of xi decouples from the equi- 
librium value for smaller x than the density of X2- This is because we have chosen a small 
value for cr 0) i and a large value for O2 . At O ,i2 = 0, Q Xl h 2 m 1.99, while fl X2 h 2 0.04. 
With increasing value of 00,12 (which parametrizes the size of the DM conversion (j2J)), Q Xl h 2 



decreases, while £l X2 h 2 increases. Around 00,12 



3.9, this order changes, i.e., fl Xl < fl X2 . At 



°"o,i2 = 4.4, we obtain the total relic abundance Q^h 2 = Q xi h 2 + £l X2 h 2 = 0.112, in accord 
with the WMAP observation tt T h 2 = 0.1126 ± 0.0036 Q. In Fig. (right) we plot Vt T h 2 as 
a function of 00,12 • We see that the DM conversion process plays an important role, as it 
has been also found in [ll|, [if], 21 , 22 



C. Three-component DM system 

As we have noticed before, the K = 3 case is possible even for a Z 2 X Z 2 symmetry if 
the decay of x% is kinematically forbidden. In this case there are nine different thermally- 
averaged cross sections, if we assume that mi > m 2 > 777.3 and m 2 + 777.3 > m \- 

<a(ii;X i X' i )v> = a 0}i x 10" 9 GeV -2 , <<r(ll;22)v>= O ,i2 x 10~ 9 GeV" 2 , 



<<t(11;33)v> = a-o,i3 x 10~ 9 GeV~ 2 , <(x(22; 33)w>= <x 0i23 x 10" 9 GeV~ 2 , 
<a(12;3X 123 )i;> = a , 123 x 10~ 9 GeV~ 2 , <cx(23; lX 231 )v>= a 0i231 x 1(T 9 GeV" 2 
<a(31;2X 312 )v> = a 0i312 x 1(T 9 GeV" 2 . 



(13) 



Eq. (|9]) then becomes 
dY 1 



dx - -0.264 [^j {<a(ll;X 1 X[)v> (YM - Y 1 Y 1 
+ <a(ll;22)v>(Y 1 Y 1 - ^J-?^) + <a(ll; 33)v>(y 1 Y 1 - 1 ^Y 1 Y 1 



Y 2 Y 2 



Y 3 Y 3 



- <«{ 12: 3A',_, )\) > - Y ?1?2 ) + <^(31;2X 312 )v>(y 1 y 3 - yYiYz 

y. 



- <a(23;ix 231 )i;>(y 3 y 2 -^y 3 y 2 



0.264 g l J 2 [^^J {<a(22;X 2 X , 2 )v> (Y 2 Y 2 - Y 2 Y 2 



<m _ 

dx 

+ <a(22; 33)v >(y 2 Y 2 - ^N^) + <<x(23; lX 231 )v >(y 2 Y 3 - y%Y 3 
+ <a{12 ] 3X l23 )v>(Y 1 Y 2 - yY x Y^j - <a(31;2X 312 )v>(Y 1 Y 3 - yYY 3 
- <(j(11;22)v>(y 1 Y 1 - J^yFi 



Y 2 Y 2 



0.264 gl' 2 1^1 {<a(33;X 3 X' 3 )v> (Y 3 Y 3 - y 3 y 3 



dY, 
dx 

+ <( r(23; lX 231 >>(V 2 y 3 - yY 2 Y^j + «r(31; 2X 312 )t; >{y x Y 3 - yYY 3 
- <cj(12;3X 123 )v>(y 1 Y 2 - ^YY 2 \ - < ( r(ll; 33>>fy 1 y 1 - ^Y l Y 1 

v y 3 j \ y 3 y 3 



<a(22;33)t;>(y 2 y 2 -ppy 2 y 2 ) I 



(14) 



(15) 



(16) 



where 1/fi = 1/mi + l/m 2 + l/m 3 . 

As a representative example we consider the following set of input values of the parame- 
ters: 



mi = 200 GeV , m 2 = 160 GeV , m 3 = 140 GeV 
cr i = 0.1 , cr 9 = 2 , cr 03 = 6 . 



(17) 



First we show the evolution of fl Xi h 2 (x) in Fig. [3] (left) for o"o,i 2 = cro,i 3 = er , 23 = ao,i 23 = 
°"o,3i2 = ^o^i = 0, which corresponds to the situation without the non-standard annihilation 
processes. Since m x > m 2 ,m 3 and the cross section a(ll; XiXi) is small in this example, 
the relic abundance of xi is large compared with that of \2 and xs- This changes if we 
switch on the non-standard annihilation processes. This is shown in Fig. [3] (right), where 
we have used o" 0j i 2 = cr 0i i 3 = o" 0j23 = 5.2, while a 0i i 23 = <7 0i3 i 2 = cr 0i23 i = to see the effects 
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FIG. 3: The relic abundance il Xl h 2 (x) (black), fl X2 h 2 (x) (blue) and Sl X3 h 2 (x) (red) as a function of 
x = fi/T = {(m^ 1 +m 2 " 1 +m 3 )T] _1 , where the input parameters are given in Eq. (jTTJ) . Left: Without the 
non-standard annihilation processes @ and ©. Right: (To, 12 = 00,13 = 00,23 = 5.2, while 00,123 = 00,312 = 
00,231 = to see the effects of XiXi XjXj type processes ©■ 

°f XiXi ^ XjXj type processes (DM conversion). As expected, the relic abundances of Xi 
and X2 decrease and drop below 0.1, while that of X3 does not change very much. 

Fig. H] shows the evolution of Q Xi h 2 (x) for 00,12 = 00,13 = 00,23 = 0, while 00,123 = 00,312 = 
°"o,23i = 5.1 to see the effects of XiXj Xk^ijk type processes (semi-annihilation). It is 
interesting to observe that the order of the relic abundances changes and Q Xl h 2 (x) first 
decreases as usual, but then starts to increase towards the freeze-out value. So, the effects 
of XiXi XjXj type and x%Xj Xk^ijk type processes are different. In the examples 
above, O ,y and ctq^ are chosen such that the total abundance Q T h 2 becomes about the 
realistic value 0.112. In Fig. [5] we show the total abundance Qj-h 2 as a function of cq ^ 
(solid) and 00,^ (dashed), where 0"o,ji parameterizes the size of the DM conversion (j2J) and 
&o,ijk parameterizes the size of the semi-annihilation ([3]) . As we can see from Fig. [5], only 
for small values of o"o,«i an d o"o,ijfe f ne effects on Qj^h 2 are different. 

Next we would like to address the question whether one can obtain a boost factor in the 
case of K = 3 within the framework of thermally produced DM, 3 where as before we assume 
that m\ > m 2 > and m 2 + > m\. Here we are interested in the situation that the 
scales that enter into the Boltzmann equations (jUj) are not extremely different, because this 
situation can easily be realized in many phenomenologically viable models. The parameter 
space for the case K = 3 is large, and therefore, to simplify the situation, we assume that 
only semi-annihilations are present. We further assume that all the semi-annihilation 
processes (|3j) produce the same SM model particle X, whose mass can be neglected. Then 
the differential energy flux $x of X produced from the DM semi-annihilations can be written 



Non-standard freeze-out history of DM in the K = 2 component DM system has been considered in 
13, 3| and a large boost factor of 0(1000) has been found. 
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FIG. 4: The relic abundance fl Xl h 2 (x) (black), il X2 h 2 (x) (blue) and il X3 h 2 (x) (red) as a function of x with 
o-o,i2 = 00,13 = o-o,23 = 0, while 00,123 = 00,312 = 00,231 = 5.1 to see the effects of XiXj XkX l]k type 
processes (|3|). 
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FIG. 5: The total relic abundance VItU 2 as a function of 00,12 (solid) and 00,123 (dashed). Except 00,12 
(DM conversion) and 00,123 (semi-annihilation) the input parameters are given in Eq. (|17p. 



as 



dE 



oc 



X 



< cr(ij; kX)v > 



(18) 



i<3 



where n i oo is the freeze-out value of rit, and Nx and Ex stand for the number and the energy 
of X in the process XiXj ~ XkX, respectively. Eq. (fl8|) suggests the following definition of 
the effective boost factor: 



BF 



x 



Si<j Y itOD Y jtOQ < crfa'; kX)v > 

Y 2 oo Y 2 ° oo <a(22;X)v>^ 



(19) 



where Y ij00 is the freeze-out value of Y{. The quantities with the superscript in the 
denominator are the reference quantities. Here we use as reference the one-component 
DM case with only the standard annihilation process X2X2 ~> XX, which means Y 2 



2,00 
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TABLE I: The effective boost factor BFx (HU) for the symmetric case, i.e., 00,123 = 00,231 = 00,312, with 
mi = 1000 GeV. 



m 2 


m 3 


CO, 123 




n X2 h 2 


n X3 h 2 


BF X 


720 


700 


12.6 


0.0433 


0.0319 


0.0372 


1.6 


940 


700 


417.0 


0.0007 


0.0007 


0.1109 


3.4 


600 


550 


42.3 


0.0431 


0.0259 


0.0439 


5.1 


840 


550 


7900 


0.0001 


0.0001 


0.1117 


5.7 



(0.11/2.82 x 10 8 )(GeV/m2). 4 If one neglects the p-wave contributions to < crv >, one finds 
cr° 2 = 2.03, 2.06, 2.09 for m 2 = 500, 700, 900 GeV, respectively, where 0° 2 is defined as 
< a(22;X)v >°= a$ j2 x 10~ 9 GeV" 2 ~ 1.16 0^2 x 10~ 26 cm 3 /s (see Eq. QH)). 

A large effective boost factor means large 0"o,yft- Since each semi-annihilation ([3]) produces 
a DM particle, the cross section for the semi-annihilation can be larger than that for the 
standard annihilation process ([T]). To get an idea on the size of < cr(ij; kX)v >, we first 
consider the symmetric case, i.e., 00,123 = 00,231 = 00,312- In Table I we give the effective 
boost factor BFx for various values of 00,123,^2 and 777.3 with the value of mi fixed at 
1000 GeV. These values are chosen such that the sum of the individual relic abundances is 
consistent with VLfh 2 = 0.1126 ± 0.0036. As we see from Table I, the effective boost factor 
BFx becomes large if the mass difference mi — 7773 is large and 7772 is close to 777.1. But BFx 
is at most of O(few), which would not be sufficient to explain the observed positron excess. 

Next we consider the asymmetric case, i.e., 00,123 7^ 00,231 7^ 00,312- The reason why we 
can not obtain a large boost factor in the symmetric case is that, as we see from Table I, 
Vj,oo x Yj,oc for i 7^ j are small. Therefore, to suppress fi X 3^ 2 an< ^ ^ ne sam e time to increase 
Q Xl2 h 2 , we consider the situation 00,123 < 00,231,00,312- In Table II we give some examples 
of the asymmetric case, which show that it is in fact possible to obtain a boost factor BF X 
of O(10). 

III. A MODEL WITH THREE DARK MATTER PARTICLES 



We extend the original radiative seesaw model of 3JJ so as to have an additional discrete 
symmetry Z' 2 . This can be done by introducing a SM singlet Majorana fermion x and a SM 
singlet real inert scalar in addition to the inert Higgs doublet 77 which is present in the 
original model. The matter content of the model with the corresponding quantum numbers 
is given in Table III. 



We use Y 2 °{ Xqc = 170) as F 2 ° c 
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TABLE II: The effective boost factor BFx (HU) for the asymmetric case, i.e., co,i23 7^ 00,231 7^ 00.312, with 
mi = 1000 GeV, m 2 = 900 GeV and m 3 = 550 GeV. 



°"0,123 


00,231 


°"0,312 


n xi h 2 


V X2 h 2 


n X3 h 2 


BF X 


48.0 


2000.0 


48.4 


0.0325 


0.0007 


0.0793 


14.5 


55.5 


65.0 


2000.0 


0.0003 


0.1118 


0.0002 


13.6 


90.0 


1000.0 


100.3 


0.0121 


0.0011 


0.0988 


13.6 


110.0 


1000.0 


125.3 


0.0091 


0.0010 


0.1022 


13.0 


110.0 


800.0 


133.0 


0.0080 


0.0012 


0.1032 


12.3 


110.0 


600.0 


145.2 


0.0067 


0.0015 


0.1039 


11.3 



TABLE III: The matter content of the model and the corresponding quantum numbers. Z2 X Z'^ is the 
unbroken discrete symmetry. The quarks are suppressed in the Table. 



field 


SU(2) L 


U{1) Y 


z 2 


^2 




2 


-1/2 


+ 


+ 


/? 


1 


1 


+ 


+ 


Nf 


1 







+ 


H= (H+,H°) 


2 


1/2 


+ 


+ 


V = (v + ,v°) 


2 


1/2 




+ 


X 


1 





+ 







1 










The Z 2 x Z' 2 invariant Yukawa couplings of the lepton sector are given by 

C Y = YZrfLil? + Y^UerjNl + Y£xK<\> + h.c , (20) 

and the Majorana mass terms of the right-handed neutrinos N% (k = 1, 2, 3) and the singlet 
fermion x are5 

£ Maj = \M k NtN c k + l -M xX 2 + h.c. (21) 

We may assume without loss of generality that the right handed neutrino mass matrix is 
diagonal and real. As far as the light neutrino masses, which are generated radiatively, are 
concerned, the last additional interaction term in Eq. (120]) has no influence. So the neutrino 
phenomenology is the same as in the original model. The most general form of the Z 2 x Z' 2 



5 A similar model is considered in 



11]. 
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invariant scalar potential can be written as 

1 



V 



m 2 H^H + m\rfrj H — m\<p 2 



+ ±X 1 (H^H) 2 + -AafoV + X 3 (H^H)(^ V ) + X 4 (H^)(^H) 
+\\,[{H^f + h.c.) + iA 6 4 + \\ 7 {rfH)<j? + \x^ V )<j> 2 
from which we obtain the masses of the inert Higgs fields: 



m 2 ± 



ml + X 3 v 2 /2 



m \ = m l + ( A 3 + A 4 + X 5 )v 2 /2 = m 2 2 + X L v 2 /2 



m 2 o = mj + (A3 + A 4 - X 5 )v 2 /2 



m 



I = m\ + X 7 v 2 /2 . 



(22) 



(23) 
(24) 
(25) 
(26) 



Here, (H) = v/y/2 is the Higgs VEV, and 77 
assumed that 



(Vr + iVi) I At this stage we have 



(H) = v/y/2, (r,) = ((f)) = 



(27) 



correspond to the absolute minimum. (The sufficient condition for the absolute minimum 
of Eq. (12"2"|) is given below.) As we can see from Table III, the cold DM candidates are 
N%, r) R , 77°, y and <f), where r) R as dark matter in the original model has been discussed in 
detail in |5JH33|. To proceed, we assume that the mass relations 



Mk » m v ±, m v o > 77^0 > m^, m x and m^o < m^ + m x (28) 

are satisfied. 6 These relations are chosen because we would like to meet the following 
requirements: 
1. \i — > e 7 



The constraint coming from /i — > ey is given by 64| 

3a 



B(fi -)• ery) 



F 2 (x) 



64ir(G F m 2 



2 ±) 2 



^2^ k Y^ h F 2 

k 



'Ml 



< 2.4 x 10 



-12 



(29) 



6(1 -x) 1 



'l-6x + 3x 2 + 2x 3 - 6x 2 In x) 



where the upper bound is taken from 65(. A similar, but slightly weaker bound for 
r — > /i(e)7 given in |65[ has to be satisfied, too. Since F 2 (x) ~ l/3x for x » 1, while 
1/12 < F 2 (x) < 1/6 for < x < 1, the constraint can be readily satisfied if M k « m v ± or 
Mfc >> m v ±. 



The possibility m^o < m r) o does not give any new feature of the model. 
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2. 

The extra contribution to the anomalous magnetic moment of the muon, = (g^ — 2)/2, 
is given by 64] 



5a, 



m, 



16n 2 m 2 ± 



'Ml 
m 2 ± 



(30) 



If we assume that | £ fc Y» k Y^ k F 2 ^ 



Ml 



M, 



£ fc TO^F 2 -£M |, then we obtain 



\5a^\ ~ 2.2 x 10~ 5 B(/x -> e7) < 3.4 x 10~ n 



(31) 



if Eq. fT29|) is satisfied, where the upper bound is taken from 66]. So, the constraint from 
a M has no significant influence. 



3. Stable and global minimum 

The DM stabilizing symmetry Z 2 remains unbroken if 
m\ < , m\ > , ml > , 

Ai , A 2 , A 6 > , A 3 + A 4 - |A 5 |, A 3 > -^(AiA 2 ) 1/2 , 



A 7 > -^(AiAe/3) 1 / 2 , A 8 > - l -(\ 2 \,/^ 



(32) 

are satisfied. These conditions are sufficient for Eq. (I27jl to correspond to the absolute 



minimum. Since m — m = A51;, a negative A5 is consistent with Eq. (128 



4. Electroweak precision 

The electroweak precision measurement requires 5jJ, |66 ] 



AT ~ 0.54 



m„± — m„o 

' 'IR 



m v ± 



m n o 



0.02 



+0.11 
-0.12 



(33) 



for mh = 115 — 127 GeV. Therefore, \m v ± — wi^oj , \m v ± — m^o| ^ 100 GeV is sufficient to 
meet the requirement. 

Then, with the assumption of the above mass relations, we look at the radiative neutrino 
mass matrix 



(M 



E 



167T 2 

ik jk 

16tt 2 



m 2 



m 2 



Ml 



In 



m 2 

Vj 



m 2 



Mi 



In 



m„o\ 2 



(34) 



m 2 /m„o \ 2 
In ' Vr 



M k 



m 2 /m ri o\ 2 
Vl In ' Vl 



M, 



M k 



for m^o , m^o < < M k 



m 2 



X5V 2 , we need EkYiiXjk £ °( 10 



-9\ 



Since (M u )ij are of order 10 -1 eV and m 2 ,,. ; it . , _ . . ,, 

for I A5 1 ^ O(0.1). Note, however, that this does not automatically imply that Y^ k \Y ik \ 2 ^ 
O(10~ 9 ), and in fact it could be much larger if we assume a specific flavor structure of Y^ k . 
If there exists another source for the neutrino mass matrix, we have to add it to 
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A. Relic abundance of dark matter 



Now we come to the relic abundance of DM. Under the assumption about the mass 
relations ( 12 8p we have to consider the following annihilation processes: 7 

•rj R rf R -H- SMs , • (f> (f> ■<->■ SMs (Standard annihilation) (35) 

*Vr Vr^ i * X X ^ 4> & (Conversion) (36) 

*Vr X ^ § v l ■, * X 4> ^ Vr, u l , • 4> Vr. ^ X v l (Semi-annihilation) (37) 




FIG. 6: Semi-annihilation diagram (left) and conversion (right). 

We have yet not specified the relative size of m x and m^. If x is lighter than <f), the 
conversion of x into <P is kinematically forbidden and the semi- annihilation in Fig. |f5] is the 
only kinematically allowed annihilation for x- So, we will consider below only the case 
m x > 1710. First we consider a benchmark set of the input parameter values: 

m v o R = 200 GeV , m x = 190 GeV , m = 180 GeV , 
m^± = m v o = 210 GeV , 
m h = 125 GeV , M 1 = M 2 = M 3 = 1000 GeV , (38) 
A 3 = -0.065 , A 7 = 0.1 , A 8 = 0.1 , \ L = -0.2 , 

E in X | 2 = 3(0.7) 2 , £ \Y^ = 9(0.01) 2 . 

k=l i,k=l 

With this choice of the parameter values we obtain 

< v(VrVri SMs)w > = 45.66 - 38.21/s , < a(<p<p] SMs)w >= 5.93 - 1.92/a; , 

< v(VrVr, H)v > = 0.46 + 0.29/x , < a{xx; <P<P)v >= + 77.18/x , (39) 

< °{xn R \ <P"l)v > = 0.02 + 0.01/x , < a{rf R (j)\ x»l)v >= 0.07 + 0.02/a; , 
<v(x4>;VrVl)v> = 0.07 + 0.04/x 



7 We neglect the coannihilations such as that of rf R with 77° and 77 ± below. 
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Y z /Y z , 

ret 

FIG. 7: Y x dependence of the relic abundances, ilrh 2 (dashed), Q v h 2 (black), VL x h 2 (blue), il^h 2 (red), 
where Y x controls the size of the semi-annihilation and conversion shown in Fig. [51 The input parameter 
values are given in Eq. (|38[) . 

in units of 10 -9 GeV -2 , and 

Q T h 2 = 0.1094 , Q^h 2 = 0.0062 , Q x h 2 = 0.0511 , Q^h 2 = 0.0521 , (40) 

where x = (1/m^o + l/m x + /T = n/T. As we see from Fig. El the size of the 

semi- annihilation and conversion is controlled by Y" fc x . In Fig. [7] we show the Y x depen- 
dence of the individual abundances, where we have varied J2k l^fc X | 2 > an d ^V^ref stands 
f° r (J2k l^fc X | 2 /3(0-7)) 1 ^ 2 . If Y x /Y* { << 1, the conversion of \ and the semi-annihilations 
X <p — y vl i X Vr ~^ v l become small, such that Q x h 2 in particular increases. 

\«) 

A 7 (Al)> 

h 

/ 

/ 

FIG. 8: Tree (left) and one-loop (right) level interactions with the quarks. 



B. Imposing constraints 

To be more realistic, we have to impose constraints from direct detection, collider exper- 
iments, and perturbativity, |Aj|, \Y£\, \Y*\ < 1, in addition to Eqs. f l29|) -f l33|) . which we shall 
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FIG. 9: The allowed regime in the Al(A7)— m^o plane for (Si = 10 , &2 = 10) GeV with m x = — GeV, 
to^ = TO^o — 20 GeV and Mk = 1000 GeV. The green and red points are for Xl and A7, respectively. 



do next. <fi and rf R have tree-level interactions to the quarks which are shown in Fig. [HI 8 In 
the following discussions we ignore the one-loop contributions such as the right diagram in 
Fig. [HI 9 The spin-independent elastic cross section off the nucleon cr{4>{ri R )) is given by 51] 

4tt U^mJ ) 



<4>{rf R )) 



(41) 



where is the nucleon mass, and / 
The cross sections have to satisfy 



a{<f>) 



0.112 



+ 



0.3 stems from the nucleonic matrix element 72 ]. 
°{rf R ) 



0.112 



< 1 



(42) 



where cxub^) is the current experimental limit on the cross section for the DM mass m. 
In the absence of x an d (ft, the lower mass region 60 GeV ^ m v o ^ 80 GeV is 



consistent with all the experimental constraints 53]. 
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But the elastic cross section 

a(rj° R ) ~ 7.9 x 10" 45 (A L /0.05) 2 (60 GeV/m„o f cm 2 with Ai £ 0.05 in this mass range 



may exceed the upper bound of the XENON100 result 73J], 11 7.0 x 10 45 cm 2 for the 



DM mass 50 GeV at 90 C.L. The higher mass region, i.e., 



m„o 



> 500 GeV with 



a{rf R ) ~ 4.6 x 10 46 (Al/0.1) 2 (500 GeV/m^) 2 cm 2 , will be significant for next generation 



XENON1T |81| or XMASS 



experiments such as SuperCDMS 

The presence of x an d <fi changes the situation. Firstly, the separation of two allowed 
regions of m v o disappears: As far as the relic abundance is concerned, m^o is allowed 



8 Direct detection of two DM particles has been discussed for instance in [g, 17, 18, 21 1. LHC signals of r\ 



dark matter have been discussed in [5lL [67 
9 There exist also one-loop corrections to rf^q 

10 There exists a marginal possibility to expend slightly this upper bound [71 1 

11 See also 0-[Z3l- 



See aJs 
rf R q U 



0,0. , 

See also [70(. 
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between 80 and 500 GeV, too, as we have seen, because x an d also contribute to the 
relic abundance. Secondly, the parameter space becomes considerably larger. To see how 
the allowed parameter space of the model without x an d changes, we consider a set of 
(Si = m^i — m v o , 5 2 = m v o — m^o ), for which the allowed parameter space without x 
and is very small. For (5\ = 10 , 5 2 = 10) GeV, for instance, there is no allowed range 
of m^o ^ 500 GeV [53]; the low mass range of m^o , for which the relic abundance Q v h 2 is 
consistent, does not satisfy the LEP constraint. Below we show how this situation changes 
in the presence of x an d 0- The LEP constraint implies that the region satisfying m^o ^ 80 
GeV and m^o < 100 GeV with 5 2 > 8 GeV is excluded Q. Therefore, for = 10 , 5 2 = 10) 
GeV we have to consider only m^o > 80 GeV. Further, to suppress the parameter space we 
assume that m x = m^o — 10 GeV, = m^o — 20 GeV and = 1000 GeV, and we scan 
m„o from 80 to 500 GeV. 

Fig. O shows the allowed area in the ~ m v% plane, where all the constraints are 

taken into account. The allowed mass range for m^o is extended as expected. The reason 
why there are no allowed points around ~ 100 GeV is the following. Since we keep 
the mass difference fixed, we have = m^o — 20 ~ 80 GeV there. This is the threshold 
regime for the process 00 — > W + W~ . So, for m^o just below 100 GeV, the annihilation 
cross section for is small because of small A7 in this range of m<^, and therefore the relic 
abundance Q^h 2 exceeds 0.12. We see that m^o = 80 GeV is allowed on the other hand. 
This allowed area exists, though A 7 is small, because around = 62 GeV the s-channel 
annihilation of becomes resonant due to = 125 GeV. For m^o just above 100 GeV, the 
annihilation cross section for is large because the channel to W + W~ is now open, so that 
VL^ti 1 can not supplement the anyhow small Q v h 2 . 

If we suppress the constraint from the direct detection, we have a prediction on the 
direct detection. Fig. [10] shows the spin-independent cross section off the nucleon versus 
the DM mass; the green area for the r\ DM and the violet area for the DM. We see that 
the the spin-independent cross sections are not only consistent with the current bound of 
XENON100 73| . but also are within the accessible range of future experiments. 



C. Indirect detection 



If DM annihilates sufficiently into SM particles, the resulting cosmic rays may be observ- 



41 



44j and in 7 



83 



86 



able. These are indirect signals of DM, and in fact excesses in e H 
have been recently reported. Indirect detection of DM has been studied within the frame- 
work of a two component DM system in 0, S, BQ, H S, Q> and also of the inert Higgs 



87H89j]. As we see from the semi-annihilation diagram in Fig. [6j the process 



model in [. 

produces only a left-handed neutrino as the SM particle. Therefore, the boost factor to en- 
hance the annihilation cross section to obtain an excess in e + is not available in this model. 
So, we are particularly interested in the neutrinos from the annihilation of captured DM in 
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FIG. 10: The spin-independent cross section off the nucleon is plotted as a function of the DM mass. The 
green and violet areas are for rj and (f> DM's, respectively, where we have used (<5i = 10 , 62 = 10) GeV with 
m x = m v o - 10 GeV, rrty = - 20 GeV and M k = 1000 GeV. 

the Sun 54 1 — 63 ] (see 0, \^ for a review, and 89] for the case of the inert Higgs model), be- 
cause (i) the semi- annihilation produces a monochromatic neutrino (E u ~ m^+m^ — m x for 
instance) in addition to those with E v ~ m rf R along with the continuum spectrum, (ii) these 
neutrinos can be observed at neutrino telescopes |9CH92| , and (hi) the evolution equations of 
the DM numbers in the Sun, which describe approaching equilibrium between the capture 
and annihilation (including conversion and semi-annihilation) rates of DM, are coupled now. 

We denote the number of DM particles 77, x and 4> in the Sun by Ni with i — rj,x an d 4>, 
respectively. Then the change of Ni with respect to time t is given by 



N v = C v - C A { m SM) Nl - C A { m o U)K ~ C A (VX MN V N X 



N v 



-C A {rj(j) xvl)N v N^ + C A (4> X Wl)N x N^ , 
C x - C A (xX <-> <t>4>)N 2 x - C A ( V x <-> (j>is L )N v N x 
+C A (v<t> xvl)N v N,p - C A {4>x <-> r]v L )N x N^ , 

- C A {<j>4> <-> SM)iVj + C A ( VV <-> #)i\^ + c A ( xx <-> H)N 2 X 
+C A (rix <-> <pv L )N v N x - C A {r](f> <-> x^l)N v N^ - C A {<\>x V^l)N x N^ 



(43) 
(44) 
(45) 



where the GYs are the capture rates in the Sun and the C A s are proportional to the anni- 
hilation cross sections times the relative DM velocity per volume in the limit v — > 0: 

2 , . v 2 



a 



4>(v) 



1.4 x io 2 VK (<) ) ( ^ 

200 GevV /a 



A 



7(L) 

0.1 



125 GeV 



-1 



rn 



0.112 



C v = 



(46) 



where the function /(m^^)) depends on the form factor of the nucleus, elemental abun- 
dance, kinematic suppression of the capture rate, etc, varying from 0(1) to 0(0.01) depend- 



18 



ing on the DM mass 60j, l61|- The annihilation rates C A can be calculated from 58 



CA» « •) = <a{ 'i;'> V> , Vij = 5.7 x 10- ( cm' , (47) 

Vij \ f^ij j 

with jJLij = 2mimj / (rrii + rrij) in the limit v — > 0. 

There are fixed points of the evolution equations, which correspond to equilibrium. Since 
at the time of birth of the Sun the numbers iVj were zero, the jVj increase with time and 
approach the fixed point values, i.e., equilibrium at which iV, assumes its maximal value. 
So, the question is whether the age of the Sun, t Q ~ 4.5 x 10 9 years, is old enough for Ni 
to reach equilibrium. We see from the evolution equations that the fixed point values are 
roughly proportional to (Cj/C ' A ) 1 ^ ! , implying that we need large capture rates Cj to obtain 
large Ni(t Q ). The annihilation, conversion and semi- annihilation rates at t — t Q are given 
by 

T{ij; .) = d l0 C A {i 3 o •)iV i (t )iV J (t ) , (48) 

where da = 1/2 and dij = 1 for i ^ j. The observation rate of the neutrinos r detect is 
proportional to T(ij; •). As a benchmark we use the same set of the input parameter values 
as ( 13 8 p . In Fig. [TT] we show the time evolution of 12 

T(SM) = C A (r]T] ^ SM)iV 2 /2 + C A {U ^ SM)^/2 , (49) 
I» = C A {r](j) v> x^N^ + C A (r] X <->• ^)N V N X + C A ( X <p r)v)N x N^ , (50) 
r(H = C A {m w)N 2 J2 , (51) 

scaled to 10 20 s -1 , as function of r = t/t@. As we see from Fig. [TTJ the rates seem to 
be saturated: r(SM) is in fact saturated, but T(u) does not reach its fixed-point value 
0.002 x 10 20 s- 1 at r = t/t Q = 1. The saturated value of T(SM) is 0.28 x 10 20 s" 1 for the 
input parameters ( |38|) . which is consistent with the upper limit ~ 2.73 x 10 21 s _1 for moM = 



250 GeV of the AMANDA-II / IceCube neutrino telescope 90j. As for the monochromatic 
neutrinos we obtain T(u) = 1.1 x 10~ 3 x 10 20 s" 1 and T(vu) = 1.3 x 10~ 7 x 10 20 s _1 . To 
estimate the detection rate Tdetect we use the formula [9i 



detect 



AP(^)r inc , (52) 



where A is the detector area facing the incident beam, P(E V ) is the probability for detection 
as a function of the neutrino energy P(E U ), and r inc = T/4itRq is the incoming neutrino 
flux, i.e., the number of neutrinos per unit area per second on the Earth (R & is the distance 



12 For the monochromatic neutrinos, i.e. we have added all the semi-annihilations, because for the 

mass values given in (1551) the neutrino energies are all closed to 200 GeV. Moreover, the first term in the 
r.h.s. of Eq. (|50[) (which counts neutrinos of + — m x — 190 GeV) is a dominant contribution with 
about 95%. 
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FIG. 11: The time evolution of the annihilation rates T(SM) and T(u), where r = t/t<r>, and the input 
parameter values are given in 



to the Sun ~ 1.5 x 10 8 km). 13 The probability P(E U ) may be approximated as the ratio of 
the effective detector length L to the mean free path of the neutrinos in the detector. For 
the neutrinos (anti-neutrinos) one finds: P{E V ^) ~ 2.0(1.0) x 10 _11 (L/km)(^( P )/GeV), 
arriving at 



r detect ~ 2.2(1.1) x 10 



-21 



L \ (E vip) \ fT_ 



A 



km / Vkm/ \^GeV J Vs 



yr 



(53) 



which implies that, for the input parameters of (|38|) . 0.05 events of monochromatic neutrinos 
with ~ 200 GeV per year may be detected at IceCube 90|, we have used: A = 1km 2 , 
L = 1km. 

0.05 events per year may be too small to be realistic. However, we would like to note that 
we have studied only one point in the whole parameter space. It will be our future program 
to implement the sophisticated method of [94] and to survey the whole parameter space. We 
also would like to note that if at least one of the fermionic DM in a multi-component DM 
system has odd parity of the discrete lepton number, a monochromatic left-handed neutrino, 
which is also odd, can be produced together with this fermionic DM in a semi-annihilation 
of DM's. 



IV. CONCLUSION 



We have considered the conversion and semi-annihilation of DM in a mult i- component DM 
system. We have found that these non-standard DM annihilation processes can influence the 
relic abundance of DM a lot, which has been recently found for two-component DM systems 
The question of whether a large boost factor can be obtained within the 



m 
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A sophisticated method to compute the observation rates at IceCube was recently developed in 94 1 
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framework of thermally produced DM in a three-component DM system has been addressed. 
It has turned out that a boost factor of at most O(10) can be obtained. 

As a concrete three-component DM system we have considered a radiative seesaw model 



of 31|, which is extended to include an extra Majorana fermion % and an extra real scaler 
boson <p. The DM stabilizing symmetry is promoted to Z 2 x Z' 2l and we have assumed that 
rf R (the CP even neutral component of the inert Higgs SU (2) j, doublet), x an d <t> are DM. 



We have shown that the previously found separation 51m3| of the allowed parameter space 
in the low and high mass regions for rf R disappears in the presence of x an d 4>- 

Finally, we have discussed neutrinos coming from the annihilations of the captured DM 
in the Sun. The evolution equations of the DM numbers in the Sun, which describe ap- 
proaching equilibrium between the capture and annihilation (including conversion and semi- 
annihilation) rates of DM, are coupled in mult i- component DM system. Due to the semi- 
annihilations of DM, monochromatic neutrinos are radiated, and the observation rates of 
neutrinos have been estimated. Observations of high energy monochromatic neutrinos from 
the Sun may indicate a mult i- component DM system. 
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